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Abstract Wc present a comprehensive set of observations of the interaction of p- 
mode oseiUations with sunspots using surface-focused seismic holography. Maps 
of travel-time shifts, relative to quiet-Sun travel times, are shown for incoming 
and outgoing p modes as well as their mean and difference. We compare results 
using phase-speed filters with results obtained with filters that isolate single 
p-mode ridges, and further divide the data into multiple temporal frequency 
bandpasses. The / mode is removed from the data. The variations of the resulting 
travel-time shifts with magnetic-field strength and with the filter parameters are 
explored. We find that spatial averages of these shifts within sunspot umbrae, 
penumbrae, and surrounding plage often show strong frequency variations at 
fixed phase speed. In addition, we find that positive values of the mean and dif- 
ference travel-time shifts appear exclusively in waves observed with phase-speed 
filters that are dominated by power in the low-frequency wing of the pi ridge. 
We assess the ratio of incoming to outgoing p-mode power using the ridge filters 
and compare surface-focused holography measurements with the results of earlier 
published p-mode scattering measurements using Fourier-Hankel decomposition. 
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1. Introduction 

The use of solar acoustic (p-mode) waves to probe the subsurface structure of 
active regions (ARs) was first proposed by Thomas, Cram, and Nye (1982). 
While efforts have been made to deduce properties of sunspots from interpre- 
tations of oscillations observed within sunspots (see reviews by Lites, 1992, 
and Bogdan, 2000), recent advances in sunspot seismology have been largely 
driven by the observations of strong influences of sunspots (and ARs in general) 
on externally impinging p modes. This includes both absorption {e.g. Braun, 
Duvall, and LaBonte, 1988; Bogdan et ai, 1993) and changes in phase (often 
characterized in terms of a change in travel time: e.g. Braun et at, 1992; Braun, 
1995; Duvall et ai, 1996). A prevalent, largely phenomenological, approach to 
exploiting the travel-time shifts (relative to travel times in the quiet Sun) to 
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model the subsurface properties of sunspots has been the characterization of the 
spot as a perturbation in the background sound speed. These types of models 
have been constructed using observations from a variety of local-helioseismic 
techniques, including Fourier-Hankel decomposition {e.g. Fan, Braun, and Chou, 
1995); time-distance {e.g. Kosovichev, 1996; Kosovichev, Duvall, and Scherrer, 
2000; Jensen et ai, 2001), ring-diagrams {e.g. Basu, Antia, and Bogart, 2004), 
and holography {e.g. Lindsey and Braun, 2005b). The discovery of travel-time 
asymmetries between waves propagating towards and away from sunspots (Du- 
vall et al., 1996) have lead to the inclusion of subsurface flows in many of these 
efforts. Travel times inferred from time-distance (TD) helioseismology have in 
particular been inverted to model flows and sound-speed perturbations using a 
variety of assumptions including Format's Principle and the ray approximation 
{e.g. Kosovichev and Duvall, 1997; Kosovichev, Duvall, and Scherrer, 2000; Zhao, 
Kosovichev, and Duvall, 2001; Hughes, Rajaguru, and Thompson, 2005), the 
Fresnel-Zone approximation {e.g. Jensen et al., 2001; Couvidat et ai, 2004), and 
the Born approximation (Couvidat, Birch, and Kosovichev, 2006). A consensus 
of many of these 3D inversions has emerged consisting of a relatively shallow 
(approximately 3 Mm deep) "slower" sound-speed perturbation above a "faster" 
sound-speed layer extending 10 Mm or more below the photosphere (see review 
by Gizon and Birch, 2005). 

These phenomcnological models have been useful as foundations for develop- 
ing both forward and inverse methods under a variety of approximations and 
assumptions (Gizon and Birch, 2005). At the same time, uncertainties about the 
degree to which the magnetic fields may contribute (in ways other than through 
associated thermal perturbations and flows) to phase or travel-time shifts, par- 
ticular in the near-surface layers, have persisted. Most local-helioseismic models 
of travcl-timc shifts, to date, do not include provisions for contributions from 
unresolved near-surface layers (near the photosphere, the typical vertical resolu- 
tion provided by observed p modes is around 1 Mm, e.g. Couvidat, Birch, and 
Kosovichev, 2006). Notable exceptions include some ID (horizontally invariant) 
structural inversions using ring-diagram analyses {e.g. Simmons and Basu, 2003; 
Basu, Antia, and Bogart, 2004). 

Some observations and inferred sound-speed models may show direct evidence 
of strong near-surface contributions to the helioseismic signatures associated 
with ARs. An early example of this is the predominantly near-surface sound- 
speed perturbation consistent with Fourier-Hankel analysis (Fan, Braun, and 
Chou, 1995). Birch, Braun, and Hanasoge (2007; Paper 2) examine the relevance 
of this particular result to a more recent modeling effort. Lindsey and Braun 
(2005b) have shown that helioseismic signatures beneath ARs obtained using 
holography largely vanish when a surface ( "showerglass" ) phase shift, empirically 
related to photospheric magnetic flux density, is removed from the data. Some 
peculiar properties of the 3D time-distance inversions have also been presented 
as evidence for surface "contamination." Korzennik (2006) demonstrated that 
an inferred subsurface sound-speed "plume" structure depends critically on the 
inclusion of observations made within a sunspot penumbra and umbra. A test of 
inversions for flows performed by masking only the umbra showed little effect of 
the mask (Zhao and Kosovichev, 2003). Ringlike regions of enhanced sound-speed 



ms.tex; 20/02/2008; 1:17; p. 2 



Holography of Sunspots: Observations 



in TD inversions of sunspots have also been examined as possible artifacts arising 
from the surface (Couvidat and Rajaguru, 2007). Surface effects in magnetic 
fields also include changes in the upper turning points {e.g. Kosovichev and 
DuvaU, 1997; Braun and Lindsey, 2000; Barnes and Cally, 2001). The observed 
reduction of p-mode amplitudes in spots has been shown to contribute to travel- 
time shifts independent of actual structural changes (Rajaguru et ai, 2006) 
as has reduced wave excitation (Hanasoge et ai, 2007; Parchevsky, Zhao, and 
Kosovichev, 2007). 

Schunkcr et al. (2005) and Schunker, Braun, and Cally (2007) found that 
travel-time shifts obtained from seismic holography in sunspot pcnumbrae vary 
with the line-of-sight angle projected onto the plane containing the magnetic 
field and the vertical direction. A similar effect has also been noted by Zhao 
and Kosovichev (2006) with time-distance measurements. A satisfactory theory 
explaining these observations remains to be constructed, but some preliminary 
suggestions include mode conversion (Schunker and Cally, 2006) or radiative 
transfer effects in combination with mode propagation asymmetries (Rajaguru 
et al., 2007). Whatever the cause, the observed line-of-sight dependence of travel- 
time shifts implies that a significant component of the shifts, at least in sunspot 
penumbrae, must be photospheric in origin. 

In ID inversions in global helioseismology {e.g. Christensen-Dalsgaard, Cough, | 
and Perez Hernandez, 1988) and ring-diagram analyses {e.g. Basu, Antia, and 
Bogart, 2004), surface effects are largely characterized by their frequency-dependcnt| 
contribution to the helioseismic signatures (mode or ridge frequencies). In con- 
trast, the observations used in 3D travel-time inversions arc typically made over 
a single, wide, frequency bandpass, and do not easily allow the assessment of pos- 
sible frequency-dependent surface terms. However, there is increasing evidence 
for frequency variations in the travel-time shifts observed in ARs (Braun and 
Lindsey, 2000; Chou, 2000; Lindsey and Braun, 2004b; Braun and Birch, 2006; 
Couvidat and Rajaguru, 2007). Braun and Birch (2006) found evidence for a 
frequency variation, at fixed phase speed, of the travel times measured in active 
regions using helioseismic holography. This variation exceeds the smaller fre- 
quency variation expected from travel-time shifts computed from a proxy sound- 
speed model, with properties similar to recent two-component 3D inversions, 
using the ray approximation. 

The observed travel-time asymmetries in sunspots (differences in travel times 
between the incoming and outgoing propagating waves) have been interpreted 
and modeled as due to flows (DuvaU et al., 1996; Zhao, Kosovichev, and DuvaU, 
2001; Zhao and Kosovichev, 2003). The shallow inflows, within the flrst 3 Mm 
below the surface, characteristic of some of these p-mode-based TD models 
appear to be inconsistent with outflows inferred from other methods including 
/-mode time distance (Cizon, DuvaU, and Larsen, 2000) and holography (Braun, 
Birch, and Lindsey, 2004). Some questions have been raised whether travel-time 
asymmetries may arise from other mechanisms, including the suppression of 
acoustic sources (Cizon and Birch, 2002; Hanasoge et al., 2007) or absorption 
(Woodard, 1997; Lindsey, Schunker, and Cally, 2007). 

Including magnetic fields in helioseismic models of sunspots appears to be a 
substantially more formidable task than constructing models that include only 
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thermal perturbations. Some progress has been made with MHD models of the 
absorption in sunspots observed from Fourier-Hankel decomposition {e.g. see 
the review by Bogdan and Braun. 1995). More recent efforts have also addressed 
the observed phase shifts {e.g. Cally, Crouch, and Braun, 2003; Crouch et ai, 
2005; Gordovskyy and Jain, 2007). It is expected that considerable advances 
in modeling helioseismic data will follow from the current development and ap- 
plication of hydrodynamic (HD) and magnetohydrodynamic (MHD) simulations 
{e.g. Jensen, Duvall, and Jacobsen, 2003; Tong et ai, 2003; Mansour et ai, 2004; 
Werne, Birch, and Julien, 2004; Benson, Stein, and Nordlund, 2006; Hanasoge 
and Duvall, 2007; Hanasoge, Duvall, and Couvidat, 2007; Khomenko and Col- 
lados, 2006; Shelyag, Erdelyi, and Thompson, 2006; 2007 Cameron, Gizon, and 
Daiffallah, 2007; Cameron, Gizon, and Duvall, 2007). 

Our primary motivation in this paper is to expand the measurements of Braun 
and Birch (2006). We hope that a comprehensive exposition of helioseismic ob- 
servations of travel-time shifts, and their dependence on p-mode properties, will 
promote and support improved modeling efforts, including the use of numerical 
simulations. As we are specifically interested in the importance of near-surface 
effects we also examine the relationship between the observed travel-time shifts 
and the photospheric magnetic field. Of particular importance is the measure- 
ment of frequency variations of the travel-time shifts, using methods similar to 
Braun and Birch (2006). However, we extend those measurements to include 
both mean travel-time shifts and travel-time asymmetries, and we determine 
spatial averages of these quantities over sunspot umbrae, penumbrae, and other 
magnetic regions. Our principle tool is surface-focused helioseismic holography 
{e.g. Braun and Lindsey, 2000; Braun and Birch, 2006), for which the travel-time 
shifts are expected to have the most sensitivity to near-surface perturbations. 
The use of surface-focus holography (described in Section 2) contrasts this work 
with other recent ( "lateral- vantage" ) holographic studies of ARs {e.g. Lindsey 
and Braun, 2005a; 2005b). In addition, to ensure a meaningful comparisons of 
our results (described in Section 3) to TD observations we use narrow annular 
pupils and corresponding phase-speed filters as discussed in Section 3.1. An 
overriding theme in our findings is a strong sensitivity of the results to the 
choice of filter and frequency bandwidth employed. To investigate this further, we 
also employ filters centered on the p-mode ridges (Section 3.2). The ridge-based 
filters allow a detailed comparison of surface-focused holography measurements 
of both travel-time shifts and absorption with published results of Fourier-Hankel 
analysis (Section 3.3). 



2. Analysis 

Helioseismic holography (HH) is a method based on the phase-coherent imaging 
of the solar interior acoustic field. It computationally extrapolates the surface 
acoustic field into the solar interior (Lindsey and Braun, 1997; 2000) to estimate 
the amplitudes of the waves propagating into and out of a focus point at a chosen 
depth and position in the solar interior. These amplitudes, called the ingression. 
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and egression, are estimated by a convolution of the surface oscilla- 

tion signal (tp) (typically the line-of-sight component of velocity observed from 
Dopplergrams) with appropriate Green's functions (Lindsey and Braun, 2000). 
For this work, the Green's functions are computed in the eikonal formulation 
(Lindsey and Braun, 1997; 2000). For surface- focused HH, the Green's functions 
represent propagators that evolve the acoustic field forward or backward in time 
from a position on the solar surface into the solar interior, and back up to the 
surface focus. To select a particular set oip modes, these functions are evaluated 
for over a chosen annular pupil. A dispersion correction, empirically determined 
from statistics obtained from measurements in the quiet Sun, is applied to the 
computation of the Green's functions (see Lindsey and Braim, 2000). 

The basis of our analysis consists of what are termed local control correlations 
(Lindsey and Braun, 2004a; 2005a). These are directly comparable to ccnter- 
annulus TD correlations {e.g. Duvall et ai, 1996; Braun, 1997). In the space- 
frequency domain, the correlation, 

C+ir) = {H+{r,,,)r{r,iy))A., (1) 

describes the egression control correlation, while 

C_(r) = (^(r,i.)i/:(r,z/))A., (2) 

describes the ingression control correlation. Here, ipi^, v) represents the temporal 
Fourier transform of the surface wave field, v is the temporal frequency, r is the 
horizontal position on the solar surface, and H-{t,v) and 7?+(r, i/) represent 
the temporal Fourier transforms of the ingression and egression respectively. 
The asterisk denotes complex conjugation, and the brackets indicate an average 
over a chosen positive frequency range Ai^. 

The primary quantities of interest are the travel-time shifts which are related 
to the phase of the correlations, 

,5r± ==arg[C±(r)]/27ri/o, (3) 

where vq is the central frequency of the bandpass Ai^. These represent travel-time 
shifts of the observed incoming (r_) or outgoing {t+) waves, as sampled by a 
chosen filter, relative to the travel times expected for the same ensemble of waves 
propagating in the solar model used to compute the Green's functions. Small 
systematic deviations of the quiet-Sun values from zero, which vary with pupil 
size and filter and are likely caused by imperfections in the Green's functions 
and dispersion correction are removed by subtracting averaged quiet-Sun values 
from the observed control correlation phases. Of interest are the mean travel-time 
shift: (5T„ican = {5t+ + (5t_)/2, and the difference (or travel-time asymmetry): 
i^Tdiff = (5r+ - (Sr_ . 

A 27 hour sequence of full disk Dopplergrams with one minute cadence, ob- 
tained from the Michelson Doppler Imager (MDI; Scherrer et ai, 1995) onboard 
the Solar and Heliospheric Observatory (SOHO), were used in this study. The 
data set starts on 1 April 2002, 21:01 UT, and includes several sunspot groups 
(NOAA groups 9885, 9886, 9887, and 9888) within a 60 x 60° Postel-projected 
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region. This area was tracked at the Carrington rotation rate and includes four 
sunspots with pcnumbral radii greater than 15 Mm as weU as other smaller 
spots. The three largest sunspots are very similar in size, with mean umbral and 
penumbral radii of 7 and 18 Mm respectively. 

The following steps summarize the general data reduction: 1) projection of the 
desired region from full-disk Dopplergrams to a Postel projection that rotates 
with a fixed Carrington rate. 2) temporal detrending by subtraction of a linear 
fit to each pixel signal in time, 3) removal of poor quality images, identified by 
a five-sigma deviation of any pixel from the linear trend , 4) Fourier transform 
of the data in time, 5) (optional) correction for the amplitude suppression in 
magnetic regions (Rajaguru et ai, 2006), 6) spatial Fourier transform of the 
data and multiplication by a chosen filter, 7) extraction of the desired frequency 
bandpass, 8) computation of Green's functions over the appropriate pupil, 9) 
computation of ingression and egression amplitudes by a 3D convolution of the 
data with the Green's functions, and 10) computation of the travel-time shift 
maps by Equations (l)-(3). 

The optional correction for amplitude suppression (step 5) involves dividing 
the amplitude of each pixel in the data by its root-mean-square value over the 
frequency bandpass. In step 6 we have used two sets of filters: phase-speed filters 
and ridge filters. Their description and the results obtained from each set are 
described in Section 3.1 and Section 3.2. For the phase-speed filters (Section 3.1) 
we compare results with and without the amplitude-suppression correction. The 
difference is relatively small and, for the ridge-filtering (Section 3.2), we use only 
uncorrected data. 

3. Observations 

3.1. Phase-Speed Filters 

The phase-speed filters used are of the type specified by Couvidat, Birch, and 
Kosovichev (2006), namely, the three-dimensional Fourier transform of the data 
(step 6) is multiplied by a function 

F(k, ly) = exp [-(27ri//|k| - w)'^/2Sw^] (4) 

where w and Sw are the mean phase speed and filter width respectively. We use 
the same set of ten filters (denoted A through J) of Braun and Birch (2006) 
with parameters listed in Table 1 of that paper. These filters are of the same 
type, but are somewhat narrower in width, than the eleven common filters often 
employed in TD analyses {e.g. Couvidat, Birch, and Kosovichev, 2006, Zhao and 
Kosovichev, 2006). Each filter is used with a corresponding pupil, over which the 
ingression and egression are evaluated. This pupil is a complete annulus defined 
so that acoustic rays at a frequency of v = 3.5 mHz reaching the inner and 
outer radii span the full width at half maximum (FWHM) of the squared filter. 
The parameter Sw is related to the FWHM by Sw = FWHM/[2(ln2)i/2]. The 
filters were chosen such that the sets of FWHM and corresponding pupils span 
a continuous range of phase-speed and radius respectively. 
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All of the phase speed filters used also remove the contribution of the / mode, 
a practice first used by Giles (2000). Our /-mode cutoff consists of a high-pass 
filter with Gaussian roll-off in temporal frequency. The position and rate of the 
roll-off varies with spatial wavenumber such that full transmission occurs at a 
frequency midway between the /-mode and pi ridges and fO^^ transmission 
occurs at the frequency of the /-mode. 

The filters are applied to. and the travel-time shift maps computed from, 
portions of the data extracted in 1-mHz-wide frequency bandpasses centered at 
2, 3, 4, and 5 mHz. We also compute travel-time shifts over a wider bandpass 
(2.5-5.5 mHz) typical of common TD measurements. The use of narrow fre- 
quency bandpasses, although frequently employed in HH, differs from typical 
applications of TD correlations. Consequently the wide bandpass measurements 
provide a useful check and basis for comparison. It should be noted, however, 
that the power spectra of solar acoustic oscillations are naturally limited in 
bandwidth. Power spectra computed using typical phase-speed filters applied to 
MDI data, without any additional temporal filtering, show a concentration of 
between 60 and 90% of the power within a 1 mHz bandwidth, depending on the 
choice of phase-speed filter. 

Some maps of travel-time shifts computed with phase-speed filters are shown 
in Figures 1 and 2. Figure 1 shows maps of (5r_ (top panels) and (5t+ (bottom 
panels), while Figure 2 shows the corresponding maps of (5T,„oan (top) and Sr^itf 
(bottom). The maps are stacked into columns of increasing phase-speed (left to 
right) and rows representing increasing frequency (bottom to top). For simplicity 
we will refer to each map by a number-letter combination denoting filter and 
frequency combination; e.g. "4B" refers to filter B applied to the frequency 
bandpass centered at 4 mHz. Because of the filter masking the / mode, and the 
decrease in acoustic wave amplitudes at frequencies below the pi ridge, only a 
subset of possible frequency-filter combinations produce meaningful correlations. 
The maps analyzed here are limited to filters 2D -2 J, 3B-3J, 4A-4J, and 5A- 
5J. 

Individual maps exhibit spatial relationships between the travel-time shifts 
and the surface magnetic flux explored below. The stacking of maps in Figures 1 
and 2 reveals several striking properties of the travel-time shifts, including fre- 
quency dependencies of St^, St-, Sr^cam and drdis at all phase-speed filters, and 
a surprising connection between the sign of the shifts and the value of central 
frequency of the filter with respect to the frequency of the pi ridge. In particular, 
positive travel-time shifts (for both incoming and outgoing waves) are observed 
exclusively in frequency bandwidths that are centered below the pi ridge, shown 
by the solid fine in Figures 1 and 2. Maps of St+, St^, (5r,„oan, and drdis for filters 
with frequencies closest to the pi ridge, i.e. 4B, 3C, and 2D (not shown), show 
positive shifts near sunspot umbra and negative travel-time shifts elsewhere in 
the ARs. At frequencies above (below) these values, the filters yield exclusively 
negative (positive) travel-time shifts. Filters E-J (not shown) exhibit trends 
similar to filter D. For these filters, negative values for i5t+ and St- are observed 
throughout the active region. Both incoming and outgoing time shifts measured 
with the larger phase-speed filters increase with increasing frequency, with the 
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Figure 1. Maps of travel-time shifts St— (top panels) and (5t+ (bottom panels) covering a 
portion of the region studied and showing sunspot group 9885. The columns of maps labeled A 
through D indicate the phase-speed filter used, while the rows indicate the frequency bandpass. 
The solid jagged line running diagonally through the panels connects the location of the pi 
ridge in the u-w domain for each filter, with the centers of the maps assigned to values of 
frequency and phase speed as indicated on the left and bottom edges of the plot. The dashed 
and dotted lines indicate the locations of the p2 and p3 ridges respectively. The map in the 
lowest-left position of the top set of panels shows a MDI continuum intensity image while the 
map in the same position in the bottom set shows a linc-of-sight magnetogram 
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Figure 2. Maps of the travel-time shifts i5rmoan (top panels) and Sr^ig covering a portion of 
the region studied and showing sunspot group 9885. The columns of maps labeled A through 
D indicate the phase-speed filter used, while the rows indicate the frequency bandpass. The 
solid jagged line running diagonally through the panels connects the location of the pi ridge 
in the u-w domain for each filter, with the centers of the maps assigned to values of frequency 
and phase speed as indicated on the left and bottom edges of the plot. The dashed and dotted 
lines indicate the locations of the p2 and p3 ridges respectively. The map in the lowest-left 
position of the top set of panels shows -Btot (see text) while the image in the same position in 
the bottom set isolates in shades of grey different portions of the AR for study (see text). 
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values of the outgoing shifts exceeding the incoming shifts. These filters (E- J) 
show AR travel-time differences ((^Tdiff ) that decrease with increasing frequency. 

As noted earlier by Braun and Birch (2006), the travel-times are non- linearly 
related to the surface magnetic-flux density. The quantity i?tot is derived from a 
MDI line-of-sight magnctogram assuming the magnetic field is the gradient of a 
potential, and used as a proxy for the total flux density. Figures 3 and 4 show 
plots of travel-time shifts against i^tot for phase-speed filters B and E respec- 
tively. For clarity, the scatter of individual pixel values is not shown (although 
see Figure 1 of Braun and Birch, 2006 for some examples). Instead, Figures 3 
and 4 show the average of the shifts derived from bins equally spaced in the 
logarithm of the flux density. Solid (dotted) lines indicate time shifts computed 
with (without) the amplitude-suppression correction discussed in Section 2. The 
effect of the correction is to decrease the travel-time shifts, especially at lower 
phase-speeds (and smaller pupils) and in regions of high flux densities typical 
of the sunspot umbrae and penumbrae, by amounts on the order of one to ten 
seconds. Vertical bars in Figures 3 and 4, and in most of the other plots shown in 
this paper, indicate the total deviation (maxima minus minima) of the averaged 
value as determined within three independent sub-regions containing the three 
largest sunspots. Thus, the bars include contributions not only from sources 
of random error, but also potentially systematic differences between sunspots. 
Typically, however, these deviations as a whole arc very small. 

The frequency variation shown for filter B (Figure 3) is typical of the results 
with smaller phase-speeds (and smaller pupils) which undergo the transition from 
positive time-shifts (at sub-pi frequencies) to negative time-shifts (at frequencies 
higher than the pi ridge) . The different frequencies shown in Figure 3a exhibit the 
three types of dependence on flux density of the mean travel-time shift described 
by Braun and Birch (2006). The travel-time differences (Figure 3b) show similar 
trends except at the highest flux densities, where the results for all frequencies 
approach zero. For filter E (Figure 4), it is observed that the frequency variations 
are in general larger for (5r_ shifts than for (5r-|_. As noted earlier, the difference 
(^Tdiff for these larger phase speeds show larger shifts at lower frequencies. This 
trend is opposite that observed with the mean shifts (Srmcan- 

As Braun and Birch (2006) note, the close relationship between (Jrmoan and 
Btot is consistent with predominately near-surface perturbations, but does not 
rule out subsurface perturbations that may very well correlate with surface flux. 
As discussed in Section 1, the variation of travel-time shifts with both temporal 
frequency and phase speed is of critical importance in understanding the depth 
variations of the underlying perturbations. In particular, Braun and Birch (2006) 
have suggested that the variation with frequency, at fixed phase speeds, of travel- 
time shifts may be a signature of surface effects (see also Paper 2). To quantify 
these variations, we compute spatial averages of the travel-time shifts over three 
types of regions characteristic of the sunspot groups. The first two types are 
sunspot umbrae and penumbrae, identified by brightness values less than 50% 
and 92% of the mean MDI continuum values respectively. The third region of 
interest (which we simply call "plage") is identified by values of i?tot above 
100 Gauss and excluding areas previously marked as umbrae or penumbrae. 
The panel in the lowest-left corner of Figure 2 illustrates in increasingly darker 
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Figure 3. Travel-time shifts for filter B: a) the mean travel-time shift, b) the travel-time 
difference, c) the incoming travel-time shift, and d) the outgoing travel-time shift as functions 
of the magnetic field Btot- Solid (dotted) lines connect averaged travel-time shifts, averaged 
over equally spaced bins in the logarithm of the flux density, computed with (without) an 
amplitude-suppression correction (see text). Red, green and blue lines indicate frequencies of 
3,4, and 5 mHz respectively. Vertical bars indicate the total deviation (maxima minus minima) 
of the averages between three independent sub-regions containing the three largest sunspots. 



shades of grey the three regions (plage, penumbrae, and umbrae) identified in 
this manner around NOAA 9885. The umbral and penumbral averages are of 
particular importance since they represent time shifts experienced by waves 
propagating through the immediate subsurface layers of a sunspot {e.g. within 
10 Mm depth below a 30 Mm diameter spot). Figures 5, 6, and 7 show the 
spatially averaged time shifts for the umbrae, penumbrae and plage respectively. 
These figures quantify many of the properties already noted in the travel-time 
shift maps; including the strong frequency variations at each fixed phase speed, 
and the changes of sign at smaller phase speeds. Figures 5-7 also include the 
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Figure 4. Travel-time shifts for filter E: a) the mean travel-time shift, b) the travel-time 
difference, c) the incoming travel-time shift, and d) the outgoing travel-time shift as functions 
of the magnetic field Btot- Solid (dotted) lines connect averaged travel-time shifts, averaged 
over equally spaced bins in the logarithm of the flux density, computed with (without) an 
amplitude-suppression correction (see text). Orange, red, green and blue lines indicate fre- 
quencies of 2,3,4, and 5 mHz respectively. Vertical bars indicate the total deviation (maxima 
minus minima) of the averages between three independent sub-regions containing the three 
largest sunspots. 



measurements made with the wide temporal bandpass (2.5-3.5 mHz). These 
values, as expected from the relative contributions of modes in the power spec- 
tra, fall largely between the results obtained in the 3 and 4 mHz bandpasses. 
Also noteworthy is that for the plage, Sr^is varies from slightly negative at low 
frequencies to slightly positive at high frequencies, and is mostly independent of 
phase speed. 

It is noteworthy that essentially all of the travel-time shifts observed in Fig- 
ures 5-7 show significant frequency variations. At smaller phase speeds, the 
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Figure 5. a) The mean travel-time shift, b) the travel-time difference, c) the incoming trav- 
cl-time shift, and d) the outgoing travel-time shift, averaged over the sunspot umbrae, as 
functions of the phase speed w. Yellow, red, green and blue lines indicate frequencies of 2,3,4, 
and 5 mHz respectively. The black dashed line indicates the use of the wide (2.5-5.5 mHz) 
frequency bandpass. Vertical bars indicate the total deviation (maxima minus minima) of the 
averages between three independent sub-regions containing the three largest sunspots. 



variations of both mean travel-time shifts and travel-time asymmetries show 
strong variations which often include a change of sign of the shifts. At higher 
phase speeds, the mean travel-time shifts also show large systematic frequency 
variations. The mean shifts observed at 5 mHz, for example, are typically twice 
the value at 3 niHz, with the difference being about 15-20 seconds. 

The effects of dispersion, including changes of ray paths as a function of tem- 
poral frequency {e.g. Barnes and Cally, 2001), may cause frequency variations 
of travel-time shifts. However, the frequency variations in the mean travel-time 
shifts observed here may be much larger than are expected for sound-speed 
perturbations inferred from recent inversions of travel-times. Using the ray ap- 
proximation, Braun and Birch (2006) computed differences of the mean travel- 
time shifts between 3 and 5 mHz, for a sound-speed perturbation similar to that 
of Kosovichev, Duvall, and Scherrer (2000), of about five to ten seconds for phase 
speeds less than 30 km s^^, and less than one to two seconds for higher phase 
speeds. The implications of these variations are best explored in the context of 
modeling (some initial efforts are addressed in Paper 2). However, it is expected 
that observations such as shown in Figures 5-7 may lead to methods for identi- 
fying and removing surface effects. There are no precedents in either global or 
local helioseismic inversions which offer any hope of including the contribution of 
unresolved near-surface structure without making use of the temporal frequency 
dependencies of the observables. 
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Figure 6. a) The mean travel-time shift, b) the travel-time difference, c) the incoming trav- 
el-time shift, and d) the outgoing travel-time shift, averaged over the sunspot penumbrae, as 
functions of the phase speed w. Yellow, red, green and blue lines indicate frequencies of 2,3,4, 
and 5 mHz respectively. The black dashed line indicates the use of the wide (2.5-5.5 mHz) 
frequency bandpass. Vertical bars indicate the total deviation (maxima minus minima) of the 
averages between three independent sub-regions containing the three largest sunspots. 



3.2. Ridge Filters 

In this section, we present results obtained using filters that isolate individual 
p-niode ridges. The use of ridge filters allows us to judge the sensitivity of the 
travel-time shifts (especially those experienced by waves near the pi ridge) to 
the choice of filter. Ridge-based filters also facilitate a more direct comparison 
with results obtained with Fourier-Hankel decomposition (see Section 3.3). Ridge 
filters have been used previously in time-distance helioseismology for /-mode 
studies {e.g. Duvall and Gizon, 2000; Gizon and Birch, 2002; Jackiewicz et ai, 
2007a; 2007b), and recently for p modes (Jackiewicz, Gizon, and Birch, 2007). 
Here we employ filters for HH that isolate the pi, p2, Pa, and pi ridges. At each 
temporal frequency, the filters have full transmission for wavenumbers spanning 
the midpoints between the desired ridge and the neighboring ridges. Sharp Gaus- 
sian roU-offs (similar to those used to remove the / mode in combination with 
the phase-speed filters in Section 3.1) remove the contributions above and below 
these wavenumbers. 

Unlike the common use of ray theory to define the radii of annuli (TD) or 
pupils (HH), there is no "standard" procedure for adopting a pupil geometry for 
ridge filters. After some trial and error, we settled on a fixed pupil for each ridge. 
We found that the results were not largely dependent on the outer pupil radius, 
but did change significantly with the choice of inner pupil radius. A choice of 
inner pupil radius smaller than roughly the horizontal p-mode wavelength of the 
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Figure 7. a) The mean travel-time shift, b) the travel-time difference, c) the incoming trav- 
el-time shift, and d) the outgoing travel-time shift, averaged over the plage, as functions of 
the phase speed w. Yellow, red, green and blue lines indicate frequencies of 2,3,4, and 5 mHz 
respectively. The black dashed line indicates the use of the wide (2.5-5.5 mHz) frequency 
bandpass. Vertical bars indicate the total deviation (maxima minus minima) of the averages 
between three independent sub-regions containing the three largest sunspots. Note that the 
vertical scale differs from those in the figures for the umbrae and ponumbrae. 



highest wavenumber present in the power spectra apparently produces undesired 
leakage of the oscillatory signal at the focus directly into the egression and 
ingression regressions. With these considerations, we adopted a set of pupils with 
radii 9-42 Mm (pi), 12-90 Mm (^2), 14-167 Mm (ps), and 17-195 Mm (^4). 
We have also experimented with varying the width of the frequency bandpass. 
The travel-time shift maps analyzed here were made with Av = 0.26 mHz, and 
were critically sampled with a frequency spacing of Ai//2. 

The results for the mean travel-time shifts, averaged over umbrae, penumbrae, 
and plage, are shown in Figure 8, while the results for the travel-time differences 
are shown in Figure 9. In light of the results using phase-speed filters (Sec- 
tion 3.1), what is most striking from Figures 8 and 9 is that both the mean and 
difference travel-time shifts using ridge filters are essentially (with a few noisy 
exceptions) always negative. This is also true along the pi ridge. Comparisons 
between the ridge-filtered and phase-filtered results are facilitated by overlaying 
the nearest phase-speed-detcrmined values on Figures 8 and 9. By "nearest" we 
mean that for a given 1-mHz frequency bandpass and radial order, the nearest 
phase-speed filter is that closest to the phase-speed of the ridge at the central 
frequency. For pi, the nearest filter combinations are 2F. 3C, 4B, and 5B. For 
P2, these combinations are 2J, 3E, 4D, and 5C. For ps, they are 31, 4E, and 
5D, and for p4, they are 4H and 5E. For p2, ps, and p4 there is very good 
agreement between phase-speed and ridge-filtered results for the mean travel- 



ms.tex; 20/02/2008; 1:17; p. 15 



D. Braun and A. Birch 



O 
CD 
CO 



to 
CD 

E 



CD 



c 
CD 

E 



20 



-20 

-40 

-60 
-80 
20 



-20 
-40 
-60 
-80 
20 


-20 
-40 

-60 
-80 
20 



-20 

-40 

-60 
-80 





n=2 



rn=3 




n=4 




2 3 4 5 

frequency (mHz) 

Figure 8. Mean travel-time shifts determined using ridge filters and averaged over umbrae 
(black lines), penumbrae (magenta lines), and plage (blue lines), as functions of frequency 
along the pi - p4 ridges. Vertical bars indicate the total deviation (maxima minus minima) 
of the averages between three independent sub-regions containing the three largest sunspots. 
The filled black and magenta circles (with 1 mHz- wide horizontal bars) represent comparisons 
of umbral and penumbral time-shifts determined from travel-time maps made with "nearby" 
phase-speed filters (see text). The circles are placed at the power- weighted average frequency 
for the given filter and frequency combination. The squares shown near 3 and 4 mHz in the 
top panel indicate umbral (black) and penumbral (magenta) mean travel-time shifts for the 
phase-speed filters (3D and 4C) which have higher values of phase-speed to those denoted by 
the filled circles (3C and 4B). 



time shifts averaged in umbrae or penumbrae, and reasonable agreement for 
the time differences. The largest discrepancies are clearly in the pi ridge, and 
(especially in the umbrae) involve a change in sign in the measurements of Tmcan 
and Tdiff between the two types of filters. 

It is noteworthy that the discrepancies between ridge-filtered and phase-speed- 
filtered results are largest for phase-speed filters with frequency bandwidths that 
are centered below pi . These are the same phase-speed and frequency combina- 
tions which produce the positive travel-time shifts seen in Figures 1 and 2. In 
contrast, adjacent phase-speed filters centered above the pi ridge apparently 
produce travel-time shifts in sunspot umbrae and penumbrae {i.e. the black and 
magenta squares in the top panels of Figures 8 and 9), which are very close to 
that observed with the pi ridge filter. Repeating the ridge-filtered measurements 
for some of the cases where these sign changes occur {e.g. near filter 3C) with 
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Figure 9. Difference travel-time shifts determined using ridge filters in the same format as 
Figure 8 



the same pupil as used with the phase-speed fiher yields yields deviations in 
both the mean shifts and travel-time asymmetry of only a few seconds (out of a 
total of 30-40 seconds) from results using the fixed pupil range stated earlier. 
Thus, the discrepancy in the sign of travel-time perturbations between the two 
types of filters is not the result of using different pupils. 

Figure 10 illustrates the extreme sensitivity of the sign of the travel-time 
shifts in sunspots near the pi ridge on the choice of filters. This figure compares 
the results of travel-time shifts computed with a commonly-used TD phase- 
speed filter (filter 1 of Couvidat, Birch, and Kosovichev, 2006 and Zhao and 
Kosovichev, 2006; hereafter TDl) with shifts computed with a pi ridge filter 
over a frequency bandpass between 3.5 and 5.5 mHz. Despite the gross similarity 
of the filtered power included in the measurements, the resulting maps of Tmean 
and Tdiff are drastically different. The travel-time shift maps made with the TDl 
filter have larger positive values {e.g. by about a factor of two in the mean travel- 
time shift and a factor of four in the travel-time asymmetry in the penumbrae) 
than results obtained with the filter 5 A shown in Figure 1, even though the 
mean phase-speed of these filters are both approximately 13 km s~^. We have 
found that incrementally increasing the width of a phase-speed filter, centered 
at 12.8 km s~^, produces maps with incrementally stronger positive travel-time 
shifts in sunspots. Based on our experience with a variety of filters, we find in 
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Figure 10. a) Azimuthally-averaged MDI power spectrum multiplied by phase-speed filter 
"1" of Couvidat, et a/. (2006). The grey scale indicates the logarithm of the power, normalized 
to the peak power. The dashed line indicates the mean phase speed of the filter {w = 12.8 
km s~^). b) The power spectrum at frequencies between 3.5 and 5.5 mHz multiplied by a pi 
ridge-filter, normalized to the peak power, c) mean travel-time shifts using the phase-speed 
filter, d) mean shifts using the ridge filter, e) travel-time differences using the phase-speed 
filter, and f) travel-time differences using the ridge filter. 



general that the requirement for producing positive travel-time shifts appears to 
be a disproportionate contribution to the correlations of wave power from the 
low-frequency wing of the pi ridge relative to the high-frequency wing. This is 
illustrated in Figure 11. This asymmetry apparently results from the fact that 
the mean phase speed of the filter (12.8 km s~^; shown by the dashed line in 
Figure 10a) falls significantly below the pi ridge. 

The reasons for the strong sensitivity, including sign-changes, of the travel- 
time shifts to details of the filter {e.g. width) and the relative weighting of the 
low-frequency wing of pi are not fully understood at this time. Thompson and 
Zharkov (2007) have presented evidence that a sign switch in mean travel-time 
shifts also apparently occurs with the application of filters centered half-way {i.e. 
in the trough) between the pi and p2 ridges. 

3.3. Comparison with Fourier-Hankel Analysis 

Here we compare the ridge-filtered (mean) travel-time shifts with the phase- 
shifts observed in sunspots using Fourier-Hankel (FH) analysis (Braun. 1995). 
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Figure 11. Cuts of the power as a function of frequency at two fixed wavenumbers illustrating 
the effect of the filters shown in Figure 10: a) power at a constant spatial wavenumber corre- 
sponding to a spherical harmonic degree of 1055. The dotted line indicates the unfiltered power, 
with the /, pi, and p2 ridges labeled. The solid line shows the power multiplied by the ridge 
filter for p\ . The dotted line indicates the power multiplied by a commonly used time-distance 
phase speed filter (TDl; see text). For comparison with the ridge-filtered power spectra, the 
phase-speed filtered power spectrum is multiplied by a factor so that the integrated power over 
the frequency bandwidth is the same as the ridge-filtered spectrum, b) power as wavenumber 
corresponding to a spherical harmonic degree of 1230. For both cases, it is clear that the 
principle difference between the two filter types is the relative weighting of the two wings of 
the pi ridge, such that the phase-speed filter enhances the contribution of the low-frequency 
wing relative to the high-frequency wing. 

To do this, wc use the pubhshcd values of phase shifts determined from Fourier- 
Hankel decomposition of waves around two sunspot groups, NOAA 5229 and 
5254, observed with Ca ii intensity images made at the geographic South Pole 
(Braun, 1995). Fortunately, the sunspots studied by Braun (1995) are similar 
in size to those included in this work. As noted in an earlier comparison with 
TD measurements (Braun, 1997), the phase shifts in FH analysis are divided by 
twice the angular frequency for comparison with travel-time shifts (a switch in 
the sign of the Hankel results is also needed due to different Fourier transform 
conventions). The results are shown in Figure 12. The agreement between the 
two sets of measurements is very respectable, despite differences in methods and 
data sets. In particular, the agreement is significantly better than an earlier 
comparison between FH phase shifts and TD travel-time shifts of Braun (1997), 
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Figure 12. The mean travel-time shift determined using ridge filters and averaged over 
umbrae (black lines) and penumbrae (magenta lines) as shown in Figure 8 as a function of 
frequency along the pi -p4 ridges. The squares represent equivalent travel-time shifts deter- 
mined by previously published FH decomposition methods applied to two sunspots, NOAA 
5229 (red) and 5254(green). 



which may be due to the poorer spatial resolution and lack of mode discrim- 
ination {i.e. lack of filtering) in that study. The agreement is especially good 
for the P3 and P4 ridges, and fair for the p2 ridge. There are clear systematic 
differences for pi, however, with the FH results indicating stronger travcl-timc 
shifts at higher frequencies than the HH results. 

Even under ideal circumstances, the two methods (FH decomposition and 
surface-focused HH) may be expected to yield systematic differences. For ex- 
ample, surface-focused HH (like TD) is primarily sensitive to the set of wave 
components that propagate to the surface at chosen locations, unlike FH de- 
composition which provides no such discrimination. At low phase speeds, for 
example, phase shifts and absorption coefficients determined from FH analysis 
include contributions from wave components passing under the sunspot and, at 
high phase speeds, include contributions from waves which refract to the surface 
multiple times in the sampling annulus. The expectation is that the results of FH 
analysis may be less sensitive to perturbations at the target than those obtained 
by using TD or HH. In light of these considerations, the agreement shown in 
Figure 12 is remarkable. 
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Maps of the ratio of egression to ingression power can be used 

in surface- focused HH to probe local emission and absorption properties. Rep- 
resentative maps of this quantity, computed from several frequency bandpasses 
using the pi (top panels) and ps (bottom panels) ridge filters, are shown in 
Figure 13. The quantity q;hh = 1 ^ l^+P/l^-P (which we denote as the HH 
absorption parameter) can be directly compared with the absorption coefficients 
determined using FH decomposition measurements, apH (Figure 14). The FH 
absorption coefficients are typically smaller than the surface-focused HH ab- 
sorption parameter in either the umbrae or penumbrae. Peak values of q;hh 
of around 0.7 are observed in sunspot umbrae as compared to typical peak 
values of apH of around 0.5. Also striking is the difference in behavior of apH 
and ckHH at high frequencies, with the absorption parameter from FH methods 
decreasing towards zero at much lower frequencies then the absorption parameter 
determined from HH. Some prior HH measurements have shown evidence for this 
behavior (Lindsey and Braun, 1999), and it has been speculated that the presence 
of surrounding emission (called "acoustic glories" ) , which may not be readily 
resolved by FH decomposition methods, may be at least partly responsible for 
the lower observed values of apH- It is highly likely that the high-frequency {e.g. 
v > 5 mHz) properties of the measured absorption parameters are determined 
not only by actual absorption mechanisms in the sunspots but by local emission 
properties as well. For example, a decrease in the local emission can in principle 
decrease the egression power, and hence the HH absorption parameter. What is 
known from measurements of p-mode lifetimes suggests that such a mechanism, 
as an explanation for all apparent absorption, is not viable for waves at lower 
frequencies, where the contribution of acoustic flux actually originating at a 
given target is expected to be only a very small fraction of the total observed 
egression (or outgoing) power (Braun, Duvall, and LaBonte, 1987). Finally, we 
urge caution regarding any interpretation of the observed fall-off of ann shown in 
Figure 14. No attempt has been made to assess (and remove) any "background" 
contributions to the measured egression and ingression {e.g. contributions due to 
locally generated oscillatory motion occurring within the pupils over which the 
ingression and egression amplitudes are evaluated). Thus, the observed decrease 
towards zero of the HH absorption measurements may be unphysical. 

4. Discussion 

It is our preference to let most of the results, as presented through Figures 1-14, 
speak for themselves. A summary of these observations would likely be either 
too lengthy or otherwise incomplete in that potentially important relationships 
not directly addressed the text may be neglected. It ought to be fairly clear from 
observations such as these, however, that the interaction between solar magnetic 
regions and acoustic waves is highly complex, and that no existing model is 
sufficient in explaining or predicting the complete range of observed behavior. In 
addition, it should be kept in mind that our observations represent measurements 
of only three active regions, with only one observable, using one spectral line, 
and along (essentially) one line-of-sight. It is known, from other observations. 
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Figure 13. Maps of the ratio of egression to ingression power determined using ridge filters 
and 0.26 niHz-wide frequency bandpasses centered at a) 2.47 mHz, b) 3.52 mHz, and c) 4.55 
mHz along the pi ridge, and d) 2.99 mHz, e) 4.04 mHz, and f) 5.08 mHz along the p3 ridge. 
At all frequencies, the areas containing sunspots (and other magnetic flux) are dark, which 
corresponds to < 1. At high frequencies (panels c and f), however, note the additional 
presence of brighter regions. 
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Figure 14. The surface-focus HH absorption parameter ohh determined using ridge filters 
and averaged over umbrae (black lines) and penumbrae (magenta lines) as shown in Figure 8 
and plage (blue lines) as a function of frequency along the pi - p4 ridges. The squares rep- 
resent values of the absorption coefficient determined by previously published Fourier-Hankel 
decomposition methods QpH applied to two sunspots, NOAA 5229 (red) and 5254 (green). 
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that dependencies of phase or travel-times shifts on these and other variables 
are now recognized, if not fully understood. Deep-focus methods (as applied to 
TD or HH analyses) are expected to provide further important constraints on 
models. We also note that wc have only briefly touched on the observations 
relevant to p-mode absorption in magnetic fields, with the expectation that we 
will return to this in further publications. 

An open question is the degree to which "surface effects," unknown or un- 
accounted for physical influences of magnetic flelds on acoustic waves, are im- 
portant in the modeling of subsurface structure of sunspots and active regions. 
As the set of observations shown here conflrm and expand upon those pre- 
sented by Braun and Birch (2006), it is worth restating the general conclusions 
derived there. Namely, the strong frequency variation of the measured travel- 
time shifts cannot be explained using standard assumptions, i.e. standard ray- 
approximation based modeling applied to sound-speed models that are typical of 
published 3D inversion results. To these unexplained frequency variations must 
now be added apparently "anomalous" positive travel-time shifts (both the mean 
and travel-time asymmetry) in sunspots, so called here because the conditions 
(i.e. p-mode properties and choice of fllter), under which they appear defy the 
expectations of standard assumptions and models. We will return to both of 
these issues, in the context of models of sound-speed perturbations, in Paper 2. 

It has been argued {e.g. Zhao and Kosovichev. 2006) that the sign-change 
of travel-time shifts with varying phase speed provides evidence for the relative 
lack of importance of surface effects for standard inferences from 3D inversions. 
Certainly, a change of sign is not a typical property of known "surface terms" 
in models of frequency shifts in structural inversions in global helioseismology 
and ring-diagram analyses. However, rather than identify the sensitivity of the 
sign of Tinoan Or Tdiff to the choice of fllter with a magnetic surface effect, these 
observations seriously raise the possibility that the positive values of these shifts 
represent an artifact, by which we mean a property that is more sensitive to 
the methods of the analysis than to actual physical conditions within or below 
sunspots. 
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